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ABSTRACT 

The accretion flow in the disk dominated state of black hole binaries has peak temperature and 
luminosity which vary together in such a way as to indicate an approximately constant emitting area. 
The association of this with the last stable orbit gives one of the few ways to estimate spin when the 
mass of the black hole is known. However, deriving this radius requires knowledge of how the disk 
spectrum is modified by radiative transfer through the vertical structure of the disk, as well as special 
and general relativistic effects on the propagation of this radiation. Here we investigate the extent 
to which differences in vertical structure change the derived disk spectra by calculating these for a 
range of different stress prescriptions. We find that at a given mass accretion rate the spectra are 
almost identical for accretion rates of L/LEdd < 0.1. The spectra are remarkably similar even up to 
the highest luminosities considered {L/LEdd ~ 0.6) as long as the stresses do not dissipate more than 
about 10 per cent of the gravitational energy above the effective photosphere. This is exceeded only by 
classic alpha disks with a > 0.1, but these models give spectral variation which is incompatible with 
existing data. Therefore, we conclude that disk spectral modelling can place interesting constraints 
on angular momentum transport, but still provide a robust estimate of the spin of the black hole. 

Subject headings: accretion, accretion disks — black hole physics — X-rays:binaries 



1. INTRODUCTION 

Black hole spin is difficult to measure as it gives a 
strong signature only very close to the event horizon. 
Currently the only way to probe this is via luminous ac- 
cretion flows, which light up the regions of dramatically 
curved spacetime close to the event horizon. However, 
to use these as a diagnostic of spin and/or to test of Ein- 
steins gravity in the strong-field limit requires knowledge 
of the velocity and emission structure of this material. 
Currently, this is best understood for an geometrically 
thin, optically thick accretion disk, where the material 
rotates in approximately Keplerian orbits (except at the 
inner boundary) and the gravitational energy thermalizes 
to produce a quasi-blackbody spectr um at each radius 
in the disk (jShakura fc Sunvaevlll973f L 

Black hole binaries (BHB) can indeed show spectra 
which look like this simple 'sum of blackbodies', though 
these are generally accompanied by a nonthermal tail 
to higher energies whose origin is not well understood. 
Nonetheless, where this tail carries only a minor frac- 
tion of the bolometric power, then the disk model gives 
a fairly good description of the spectral shape. More 
compellingly, a sequence of such disk dominated spec- 
tra from the same object at different mass accretion 
rates shows that the bolometric disk luminosity (span- 
ning factors of 10-5 0) relates to the ma ximum tempera- 
ture as Lboi oc T^ ax (jKubota et al.l200"ll ). or equivalently, 
that the disk has an approx imately constant inner radius 
(jEbisawa et al.lfl99lL Il993f ). This observation is exactly 
the behavior predicted by General Relativity, as the last 
stable orbit forms a fixed size scale for a geometrically 

1 Department of Physics, University of Durham, South Road, 
Durham, DH1 3LE, UK 

2 School of Natural Sciences, Institute for Advanced Study, Ein- 
stein Drive, Princeton, NJ 08540 

3 Chandra Fellow 



thin disk (see e.g. the review bv lDone et al.ll2007L here- 
after DGK07). 

However, these disk dominated spectra are surprising 
in the context of the Shakura-Sunyaev disk models as 
they are seen at high luminosities, typically 0.1 — 0.5 
-^Edd (e.g. DGK07). The classic prescription in which 
the viscous stress is directly proportional to total (gas 
plus radiation, P to t = -Pgas + -Prad) pressure, is thermally 
and viscously unstable where radiation pressure domi- 
nate s, i.e. for luminosities above ~ 0.05LEd d for BHB 
(e.g. iLightmar] 11974 ; IShakura fc Sunvaevlll976h . Above 
this luminosity, the inner regions of alpha disks exhibit 
limit cycles, in direct confl ict with the observed stabil- 
ity of the disk spectra fe.g. lNowaklll995l ). Yet the alpha 
prescription does seem to be able to reproduce the long 
timescale outburst /quiescence behavior of the outer, gas 
pressure dominated re gion in the d isk instability model 
(see e.g. the review bv lLasota|[200l . 

Thus it seems likely that the stress is approximately 
proportional to gas pressure where gas pressure dom- 
inates, but that it scales less sensitively with temper- 
ature than predicted by radiation pressure when radi- 
ation pressure dominates. Ultimately, the question of 
the stress scaling will be answered by numerical sim- 
ulations of the Magne to-Rotational Instability (MRI, 
iBalbus &: Hawlevlll991h which forms the physical basis 
for "effective viscosity". However, this is beyond the 
scope of current codes as it requires coupled radiative- 
magneto-hydrodynamic simulations in full general rela- 
tivity, covering a large radial extent but simultaneously 
resolving the small scale height of the thin disk!. Un- 
til these become available, ad hoc stress scalings are the 
only way to model the structure of a radiative accretion 
disk. 

The choice of stress prescription can affect the spec- 
trum as well as the stability of the disk as it determines 
the surface density, and therefore, the overall optical 
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depth of the disk. This emission thermalizes to a true 
blackbody only if the disk is effectively optically thick to 
absorption at all frequencies. Free-free (continuum) and 
bound-free (photo-electric edge) absorption opacity both 
drop as a function of frequency while electron scatter- 
ing is constant, so the highest energy photons from each 
radii are less likely to thermalize. The spectrum then 
becomes a modified blackbody whose frequency depen- 
dence can, in general, differ significantly from Planckian. 
However, the effects of Compton scattering tend to yield 
a Wien-like tail and the hottest parts of the disk are 
generally well modelled by a color-corrected (or diluted) 
blackbody, with effective temperature which is a factor 
/col (termed a color-temperature correction) higher than 
for complete thermalization. The full-disk spectrum is 
then a sum of modified black bodies, but this can likewise 
be approximately described by a c olor-temperature cor- 
recte d blackbody disk spectrum (|Shimura fe Takaharal 
Il995h . This color-temperature correction depends on the 
vertical temperature and density structure of the disk, 
both of which can only be fully determined once the 
stress is specified. 

Since the observations show to zeroth order that L oc 
T 4 , the color-temperature correction must stay approxi- 
mately constant (the alternative, that it changes in just 
such a way as to mask any change in innermost extent of 
the disk, violates Occam's razor). This used to be con- 
troversial, with different models of al pha disks (with only 
H and He) giving different results. IShimura fe Takaharal 
([l99lih found that / co] remains rel atively consta n t over 
a wide range in luminosity, while iMerloni et al.l (2000) 
had constant / CQ i only at > O.liEdd, increasing be- 
low (jGierliriski fe Donell200l hereafter GD04). Since 
it is physically rather unlikely that absor ption becomes 
less e ffective at lower temperatures, the IMerloni et all 
(2000) result seems rather to be an artifact of the as- 
sumed constant vertical density structure. However, 
both these calcula t ions a re now superseded by the mod- 
els of iDavis et al.l (|2005l ). which includes fully non-LTE 
calculations for all abundant elements together with 
the self-consistent radiative transfer and disk structure. 
These indeed show that / co i is approximately constant for 
an alpha disk, but only below temperatures of ~ 1 keV. 
Above this, the disk becomes so ionized that even photo- 
electric absorption opacity becomes negligi ble, so / co i in- 
creases markedly (|Davis et al.l 12005. 2006). A constant 
/col can be recovered even at these high temperatures by 
simply decreasing a to ~ 0.01, but this is inconsistent 
with the a ~ 0.1 required in order to matc h the rapid 
rise t o outburst seen in transient systems (e.g. lKing et al.l 
12001 . and still gives a disk which should undergo limit 
cycle oscillations when radiation pressure becomes dom- 
inant. 

Thus it seems that both spectra and stability of the 
observed disk spectra in BHBs are inconsistent with 
an alp ha stress pr e scri ption. Here we use the meth- 
ods of IDavis et all (|2005D and IDavis fe Hubenvl ()2006l ) 
to examine two alternative models, where the stress 
is directly proportional to gas pressure alone, here- 
after called the beta disk, or to the geometric mean of 
the g as and total pressure, hereafter called the mean 
disk. (ISakimoto fe Coronitilll98lUStella fe Rosnerlll984l : 
ITaam fe Linl 11984 lLaor fe Netzerl Il989t iMerlonil 120031 ). 
When we fix a = 0.1 to agree with requirements of the 



disk instability model, we find that either of these give 
a better match to the spectral constraints than an alpha 
prescription as the resulting disk is denser and more op- 
tically thick. Thus these alternative stress scalings give 
disks which are able to maintain a relatively constant / co i 
even at high temperatures, providing a better match to 
the observed spectra as well as the stability properties. 

While this gives important insight into the approxi- 
mate form of the MRI stresses, we show that all stress 
prescriptions considered here (alpha, beta, and mean) 
give the same f co \ below ~ 1 keV, suggesting that disk 
spectra can be used to give a robust estimate of the inner 
radius of the disk, hence observationally constraining the 
black hole spin. 

2. ALTERNATIVE STRESS PRESCRIPTIONS 

The m agnetic nature of the viscosity has long been rec- 
ognized (IShakura fe Sunvaevlll973t lEardlev fe Lightmanl 
119751 : Ilchimarulll977l : iGaleev et al.lll979T). motivating al- 
terna tive stress prescriptions (e.g. iSakimoto fe Coronitil 
1981). The two most popular are where stress is di- 
rectly proportional to gas pressure alone (beta disks, 
iStella fe Rosnerl 119841 ) or to the geometric mean of 
the gas and total pressur e , A / P gas -Ptot (mean disks, 
ITaam fe Lirl Il984t IMerlonil [2003h . Beta disks are sta- 
ble at all luminosities below the Eddington limit (above 
which the equations break down as winds and/or advec- 
tion become important), while the mean disks show limit 
cycle s above ~ 0.3LEdd for stellar remnant black hole s 
(e.g. iHonma et ail fi99ll ; IMerloni fe Navakshinl [2006h . 
The mean disk is thus inconsistent with the observation 
that stable disk spectra exist up to at least ~ 0.7LEdd 
(DKG07) , but the onset of the instability is very sensitive 
to the exact scaling of the stress. For the more general 
prescription which is proportional to -P|Ls-Ptot ^ then the 
obser ved stability limit requires p ~ 0.56 (Kat o et al.l 
1998), which is negligibly different to the mean disk 
(p = 0.5) in terms of steady state disk structur e. 

We use the bhspec code of IDavis fe Hubenvl (|2006l ) to 
calculate the self-consistent vertical structure and result- 
ing spectra for three stress prescriptions: alpha, beta, 
and mean disks. We consider a geometrically thin disk 
and break it up into a series of logarithmically spaced 
annuli. We then calculate the one-dimensional vertical 
structure and radiative transfer in eac h annulus using 
the TLUSTY stellar atmosphere code (|Hubenv fe Land 
1995). The total spectrum is then integrated over the 
disk surface, accounting for the effects of space-time cur- 
vature on the photon geodesies. Three parameters de- 
termine the structure of each individual annulus: the 
radiative flux at the disk surface; the surface density, E; 
and a parameter which determines the strength of the 
tidal gravity. 

Since mechanical energy dissipates directly within the 
disk "atmosphere" , the radiative flux F varies with 
height z in the disk, requiring an additional assump- 
tion to uniquely determi ne the structure. Following 
IShakura fe Sunvaevl (j!973l ). we assume that the dissipa- 
tion dF/dz is locally proportional to the density p. Using 
the definition of column mass, dm = —pdz, this is equiv- 
alent to 

dF__2aT^ 

dm S ' W 
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Fig. 1. — Luminosity-Temperature relations derived from fits of bhspec with diskbb as described in J3] The left panel corresponds to fits 
using a CCD (Swift XRT) response while the right panel utilizes a proportional counter (RXTE PCA) response. We consider four different 
stress prescriptions: alpha disks with a = 0.01 (black) and a = 0.1 (red), beta disks with a = 0.1 (blue), and mean disks with a = 0.1 
(green). These all overlap apart from the alpha disks with a = 0.1 at high luminosity, showing how robust the spectra are to different 
stress prescriptions as long as the disk remains optically thick. The solid black lines show the L — T relation expected from black hole with 
a* = 0.5 for a color-temperature correction of 1.6, 1.8, 2.0 and 2.2 from top to bottom, respectively. 



where T e ff is the effective temperature of the annu- 
lus. Note that the right hand side terms, and there- 
fore dF/dm are independent of height with this assump- 
tion. We cannot use the stress prescription to specify 
a local viscous-like dissipation as this is proportional to 
some combination of radiation and gas pressure. These 
pressures are largest at the disk midplane due to hydro- 
static equilibrium, so the dissipation is also largest there, 
and the resulting temperatures are so high as to lead to 
the disk becoming Rayleigh- Taylor unstable. This would 
lead to convection, presumably driving the disk struc- 
ture towards the marginally stable condition where the 
dissipation is proportional to density. Time and horizon- 
tal averages of sh earing box simulations ([Turner! 12004 
iHirose eta l. 2006) indicate that a somewhat greater frac- 
tion of the dissipation occurs at low column near the 
surface. 

We use a vertically integrated disk model to determine 
the above parameters for each annulus. The tidal gravity 
is specified by the choice of spacetime and radius, while 
the radiative flux at the surface is i ndependent of stress 
presc ription for a thin disk (see e.g. IShakura fc Sunvaevl 
1973). Therefore, we only require a stress prescription 
to determine the surface density in this model. (Note, 
however, that the theoretical arguments which motivate 
the stress prescription may also require modifications to 
the dissipation profile, a point we will discuss further in 
f|5]) The surface density may be determined by enforcing 



angular momentum conservation 



r r(f> dz = 



MnD{r) 
~2^~ A(r) 



(2) 



where t t $ is the accretion stress, M is the accretion rate, 
and O is the Keplerian frequency. The functions A{r) 
and D(r) (both — > 1 fo r radii much larger tha n the last 
stable orbit) defined in Rif fert fe Heroldl (|1995l ) incorpo- 
rate the no-torque (stress free) inner boundary condition 
and parameterize the deviation of general relativistic disk 
structure from the Newtonian limit. 

The integral in eq. {2} presents a computational dif- 
ficulty since r r ^ depends on local variable such as tem- 
perature and density for the stress prescriptions consid- 
ered here. Therefore, the exact computation requires the 
full vertical structure which, in turn, requires a choice 
of S. Therefore, an exact solution can only be com- 
put ed via iteration p rovided an initial guess for E (see 
e.g. lDavis et aLlfeOOSh . Given the approximations inher- 
ent in our stress prescriptions, we instead solve for S us- 
ing moderately crude vertical averages to replace J r r ^dz 
with 2H{T rr j,). Here H is an approximate vertical scale 
height for the disk determined via the equation of hydro- 
static equilibrium 



= pn^W 



(3) 



dz r "~B(r)' 

where B(r) and C(r) are relativistic factors (again both 
— > 1 fo r radii much larger than the last stable orbit) de- 
fined in lRiffert fc Heroldl (| 19950 . This is solved by substi- 
tuting z — > H and d/dz — > H -1 , We then specify (P ra d) 
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Fig. 2. — The upper, middle and lower panels show the effective optical depth to the disk midplane, total surface density and effective 
temperature as a function of radius, respectively. We compare alpha disks with a = 0.01 (black, dot-dashed) and a = 0.1 (red, solid), beta 
disks with a = 0.1 (blue, dotted), and mean disks with a = 0.1 (green, dashed). T e fj is independent of stress prescription, (a) shows these 
disk parameters for logZ = —1 while (b) has log / = —0.2. 



and (Pgas) in terms of E, yielding two equations EJ) f° r 
two unknowns (S, H). The resulting E evolves smoothly 
between the gas and radiation pressur e dominated lim- 
its. Previous work (jDavis et al.l l2005h has shown that 
the resulting estimates for E typically agree with a more 
precise iterative method to < 30%. 

Our results also agree with the gas and radiation pres- 
sure dominated limits which h ave been worked out in de- 
tail by previous authors (e.g. IShakura fc Sunyaevill973l; 
Stella fe Rosnen Il984t lLaor fe Netzerl Il989t iMerlonil 
2003T ). We give the asymptotic behavior at large radii 
(away from the inner boundary condition) below as this 
is useful in understanding the subsequent plots. For the 
gas pressure dominated limit all three stress prescriptions 
give identical results where 

E gas oc m^ 5 l^ 5 r~ 3 ^a~ 4/5 (4) 

where I = L/LEdd is the Eddington ratio, r = R/R g 
is radius in unit of the gravitational radius, and m = 
M/Mq is the black hole mass. For the radiation pressure 
dominated limit the surface density is sensitive to the 
stress prescription, and for our choices gives 

Erad.aCxZ-V 3 /^- 1 , 

S^ocm 1 / 5 * 3 / 5 ^ 3 / 5 ^ 4 / 5 , 
S rad , ncan ocm 1 /9r 1 /9 r V3 a -8/9^ (5) 

Since the BHBs considered here cover only a small 
range in mass, the weak dependence of E on m can be 
ignored. In what follows, the Eddington ratio depen- 
dence will be most important. For the alpha and mean 
disks E decreases as I increases, but with a much weaker 



dependence for the mean disk. In contrast, E increases 
with I in the beta disk model. The susceptibility of the 
stress prescriptions to viscous instability follows directly 
from these relations. The condition for stability is that 
dE/dM > (jLightman fc Eardlevj|l97l . Since M oc I, 
we see that only the beta disk is truly stable in the limit 
where radiation pressure dominates entirely, though the 
mean disk is also very close to the stability condition. 

3. LUMINOSITY- TEMPERATURE PLOTS AND THE 
COLOR-TEMPERATURE CORRECTION 

We now use this bhspec code (jDavis fc Hubenv![2uu6L 
see also 5j2| to investigate the spectra produced by these 
different stress prescriptions for a black hole of 10 M©, in- 
clined at 60°, with spin of a* = 0.5. We use the XSPEC 
spectral fitting package to convolve the models with the 
detector response from current X-ray instruments in or- 
der to simulate what would be observed. We do this 
for accretion rates equally spaced in log? from -1.4 to 
-0.2 through both the RXTE PCA (3-20 keV) and Swift 
XRT (0.1-10 keV) responses. These two instruments are 
broadly characteristic of any proportional counter (e.g. 
GINGA, ASTROSAT) and CCD (e.g. Chandra, XMM- 
Newton, Suzaku XIS) responses, respectively. We fix the 
source distance at 8.5 kpc (Galactic center), and assume 
a Galactic absorbing column of Nh = 1.2 x 10 21 cm -2 
(approximately the lowest possible at this distance) . We 
note that for log/ = —0.6, this gives a very good esti- 
mate of the disk dominated spectrum of XTE J1817-330 
at the peak of its outburst (spectrum 001: iRvkoff eTall 
[20071 : iGierlinski et al.ll2008l hereafter GDP08). 

We also include a small Comptonized tail, with bolo- 
metric luminosity of ~ 5% of that of the disk, with pho- 
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Fig. 3. — Color-temperature corrections from fits of bhspec with kerrbb as described in [|3] The left panel corresponds to fits using 
a CCD (Swift XRT) response while the right panel is for a proportional counter (RXTE PCA) response. We consider four different 
stress prescriptions: alpha disks with a = 0.01 (black) and a = 0.1 (red), beta disks with a = 0.1 (blue), and mean disks with a = 0.1 
(green).. Comparison with Fig 1 shows that the results for the color-temperature correction are not the same as for diskbb, where the CCD 
bandpass gives an almost constant / 00 ; = 1.65 — 1.75 as opposed to the change from 1.6-1.9 seen here. The proportional counter likewise 
has f co i = 1.7 — 2.0 with diskbb compared to 1.4 — 1.8 here. 



ton index of T = 2.2 and electron temperature fixed 
at 100 keV (i.e. outside the bandpass of either instru- 
ment). The seed photons are assumed to be a blackbody 
at the maximum inner disk temperature, set by noting 
that disk fits to the peak spectrum of XTE J1817-330 
give kT disk ~ 0.9 keV (GDP08). We use this as our fidu- 
cial seed temperature for log / = —0.6, and then change 
it by a factor (10 01 ) 1 / 4 = 1.059 for each 0.1 dex change 
in logZ. We also add 1% systematic uncertainties to our 
simulated PCA spectra as is usual with real data. 

These simulated spectra are then fit wit h a very simple 
multi color disk blackbody model, diskbb (Mi tsuda et al.l 
1984). This assumes that the disk temperature T(r) oc 
r~ 3 / 4 , i.e. has continuous stress at the inner boundary. 
The emitting area can still be deriv ed from the data by 
applying a correction factor for this (|Kubota et al.|[l998l) 
as the overall shape of this spectrum is very similar to 
that produced by a str ess-free inner boundary condition 
(|Gierlmski et alJ[200l . We also include a Comptonized 
tail produced from (blackbody) seed photons tied to the 
inner disk temperature and absorption fixed at Galactic 
values. The disk so dominates the CCD bandpass that 
the Compton tail is poorly constrained, and allowing the 
spectral index to be free gives best fit parameters, for 
which the Comptonized tail is very steep. This is due to 
the bhspec disk shape being subtly different to a simple 
sum of blackbodies, firstly as relativistic effects smear 
out the Wien tail and secondly as the color-temperature 
correction is not constant as a function of radius. These 
combine to give a broader disk spectrum than predicted 



by diskbb so the low temperature Comptonization is re- 
quired to extend the emission to higher energies. Instead, 
we make the assumption that simultaneous high energy 
data exist to constrain the spectral index (such as from 
an RXTE monitoring campaign: e.g. GDP08), hence we 
fix this at the input value of 2.2 in the CCD fits, but 
allow it to be free in the proportional counter models. 

Fig. [T] shows the resulting L — T diagram from fitting 
the CCD and proportional counter simulations for the 
two new stress prescriptions considered here (green: beta 
disk with a — 0.1, blue: mean disk with a — 0.1), to- 
gether with the two previous alpha disk models for com- 
parison (black: a — 0.01 and red, a — 0.1). All the stress 
prescriptions which yield large surface density (a = 0.01, 
mean and beta disks) show very similar behavior. Only 
the alpha disk with a = 0.1 gives significantly higher 
temperatures at high luminosities. However, at low lumi- 
nosities, all the stress prescriptions give a fairly g o od T 4 
relation, showing clearly that the iMerloni et alj (2000) 
result of a substantial increase in color temperature at 
low luminosities is indeed an artifact, and is not repre- 
sentative of a gas pressure dominated beta disk (c.f. the 
discussion in iGierliriski fc Pond 120041 where the Mer- 
loni et al. result was used to argue against a beta disk 
model). 

The luminosity in these plots is derived as in the L — T 
plots of GD04 and DGK07 by correcting the disk flux 
and temperature for inclination angle and relativistic ef- 
fects. The general relativistic corrections are taken from 
Zhan g~et al.l ([1997) for a Schwarzschild black hole as this 
is closer to an a* = 0.5 disk than the alternative tabula- 
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Fig. 4. — a) The best-fit diskbb models to simulated bhspec disk spectra with a* = 0.5 and logZ = —0.3 for the proportional counter 
(red, solid curve) and CCD (green, dashed curve) responses. For comparison, we also plot actual bhspec model (blue, dotted curve) used to 
generate the simulated spectrum. The diskbb spectrum is narrower than bhspec so the best fit models adjust as the bandpass changes. The 
CCD gives a lower temperature by a factor 1.15, and the discrepancy at 5-10 keV is taken up by the freedom in fitting the Comptonized 
tail, b) A similar comparison of kerrbb models with the same bhspec model plotted in panel a. However, in this case the kerrbb spectra 
are not best fit models. Instead, they have spin and mass accretion rates set at the intrinsic parameters for the bhspec model, but with 
color-temperature corrections chosen to have the same ratio (1.15, see text for further details) as the best-fit diskbb temperatures. It is 
plain that the color-temperature corrections derived from the diskbb fits are not appropriate for kerrbb because the models have different 
spectral shapes. 



tion for a» — 0.998. However, these have very little effect 
at 60°, as this is where the Doppler blueshift and gravi- 
tational redshifts approximately cancel, thus the flux to 
luminosity conversion is very close to a simple disk area 
correction of L — 2ttD 2 F/ cosi = 4ttD 2 F for cosi = 0.5, 
and the temperature correction is negligible. However, 
the recovered flux is less than the input value (e.g. for 
the highest luminosity points made from log/ = —0.2, 
while the dense disks give —0.26) as limb darkening is 
present in the bhspec radiation transfer. 

The solid lines show the predicted L — T relation for a 
constant color temperature of 1.6, 1.8, 2.0 and 2.2 again 
as in GD04, but scaling the disk area to that expected 
from an a* = 0.5 black hole. Clearly the CCD data are 
very close to a constant value of / co i = 1.7, while the 
proportional counter bandpass gives a small change in 
/coi from 1.8-2.0 for the dense disks, and 1.8-2.2 for the 
alpha disk with a = 0.1. 

The top panels of Fig. [2] show the effective optical 
depth of the disk (r off = \/r a bs(T cs + r abs )) as a func- 
tion of radius for log/ = — 1 and —0.2 respectively. For 
comparison we also plot the surface density (middle) and 
effective temperature (bottom) for the same models. For 
alpha disks the minimum in r e s closely corresponds to 
the minimum in E. Furthermore, these minima are close 
to the radius of maximum T g, so the majority of the flux 
will be produced in the regions of the disk which have the 
lowest opacity and hence the highest color temperature. 
Comparison of the top and middle panels show that T e g 
is very similar in shape to E. For a given alpha, the beta 
disk always yields a larger E (and, therefore, a larger 
r ff) than the alpha disk, with the mean disk midway 



between them. Due to the different scalings with / (see 
Eqs. [4] and [5]), the models become more discrepant as 
luminosity increases. The alpha disk with a = 0.1 (red, 
solid curve in Fig. [2|) always has the lowest effective op- 
tical depth, and so gives the largest color-temperature 
correction. However, the effect of this is small at low 
luminosities (Fig. [T]) as there is still substantial absorp- 
tion opacity (Fig. [2^) and the heating mainly occurs 
at high effective optical depth. Thus, the photosphere 
acts simply as an atmosphere. By contrast, at the high- 
est luminosities, the alpha disk with a = 0.1 approaches 
r c ff = 1. All of the flux is dissipated above the effective 
photosphere itself, and the color-temperature correction 
becomes much larger. 

We repeat the spectral fitting, replacing the simple 
diskbb model by the more physical kerrbb spectrum 
(|Li et al.l 12005). Like diskbb, kerrbb is a multicolor 
blackbody model, but it includes all of the relativistic 
effects modeled in bhspec i.e. all the special and general 
relativistic smearing of the intrinsic disk emission due to 
rapid rotation in strong gravity. These broaden the spec- 
trum so that it is not so sharply peaked as predicted in 
the diskbb models, but unlike bhspec it assumes that 
the intrinsic emission from each radius is just a color- 
temperature corrected blackbody, and that this color cor- 
rection factor is the same at all radii. Therefore, the 
derived / co i will arise solely from diff erences in the treat - 
ment of the surface emission (see e.g. lShafee et~aTl l2006). 
We fix the mass, distance, inclination and spin, and fit 
for the mass accretion rate and color-temperature cor- 
rection assuming no stress on the inner boundary and no 
returning radiation, but including limb darkening. Fig. 
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[3] shows these for the CCD (left panel) and proportional 
counter (right panel) data respectively. Again the three 
dense disk prescriptions all show similar results, with the 
alpha disk with a = 0.1 giving higher color-temperature 
correction at high luminosity, as before. However, the 
values of the derived color-temperature corrections are 
different to those derived from the diskbb approach. 
The color-temperature correction for kerrbb does not 
stay constant in the CCD bandpass, apparently in con- 
flict with the observed L oc T 4 relation from diskbb. It 
increases from 1.6—1.9, which at least goes through the 
/col = 1-7 constant value seen with diskbb. By contrast, 
the proportional counter bandpass gives values of the 
kerrbb color-temperature correction which are always 
lower than those seen from diskbb, as well as spanning 
a wider range, from 1.35-1.8 

The answer to these issues lies in the detailed shape 
of the spectrum. Fig. [4^ shows the initial bhspec 
disk (blue, dotted line) and its best fit diskbb spectra 
to the proportional counter (red, solid line) and CCD 
(green, dashed line) bandpasses for one of the dense 
disk prescriptions (a = 0.01) for an input luminosity 
of log/ = —0.3. The higher energy bandpass of the 
proportional counter weights the fit to higher temper- 
atures/lower normalization than those for the CCD, and 
the diskbb model is a good fit to the bhspec above the 
3 keV lower limit of the proportional counter bandpass. 
However, the best fit to the CCD data does not look so 
compelling. It is clear that the diskbb spectrum is nar- 
rower than the bhspec disk emission, and the 'best fit' 
matches this only in the 1-3 keV range where the signal- 
to-noise is maximum. At higher energies, the Compton 
tail provides additional freedom to match the data. 

Fig. 2)3 shows the kerrbb models with mass accretion 
rate chosen to give log I = —0.3 for a black hole with a* = 
0.5 (i.e. 8.537 x 10 18 g/s), with color-temperature correc- 
tion of 1.7 for the proportional counter (red). The ratio 
of best fit temperatures from the diskbb models above 
implies that the color-temperature correction should be 
a factor 1.15 smaller in the CCD compared to the pro- 
portional counter i.e. giving f co i = 1.47 for the CCD 
(blue). However, the kerrbb spectrum is much broader 
than diskbb as it includes the relativistic smearing of the 
continuum, and the proportional counter model is actu- 
ally a fairly good description of the bhspec disk shape 
apart from around the peak at ~ 1 keV. This is because 
bhspec includes the photoelectric edge features in the at- 
mosphere, calculating the radiative transfer through the 
ion populations rather than simply assuming a diluted 
blackbody. These atomic features, from partially ion- 
ized oxygen/silicon/iron L shell, depress the continuum 
around 1 keV from that predicted by the kerrbb models. 

Thus the kerrbb model appropriate for the propor- 
tional counter data gives a good match below ~ 0.6 keV 
but overestimates the flux at the peak. Decreasing the 
color-temperature correction simply increases the dis- 
parity between the kerrbb and bhspec disk spectra in 
the CCD bandpass, especially in the crucial 1-3 keV re- 
gion where the signal-to-noise is high. Thus the appro- 
priate kerrbb model would have slightly higher color- 
temperature correction for the CCD than for the pro- 
portional counter, as observed (see Figs. 3a and b). 
These differences in spectral shape mean that color- 
temperature corrections derived from diskbb cannot be 



simply applied to kerrbb (and vice versa). 

4. COLOR-TEMPERATURE CORRECTION AND 
EFFECTIVE OPTICAL DEPTH FOR ALL SPIN 

We use the insights derived above from the particular 
case of a* = 0.5 to extend our analysis of r e ff and / co i to 
a larger range of black hole spins and luminosities. We 
are primarily interested the role T e g plays in determining 
/coi. Since our fitting procedure is most sensitive to the 
hottest annuli, we will focus only on T e g at the radius 
where T c g is maximum. For each a* and / we calculate 
the radius of maximum T e g , and evaluate T e g at the mid- 
plane for the four stress prescriptions considered in 
Fig. [5] shows contour plots of the effective optical depth 
at the radius of maximum flux (c.f. Fig. [2]) as a function 
of spin and luminosity for four different stress prescrip- 
tions. The calculations span a range of a* from to 0.8 
and log / from -2 to 0. (We consider higher spins below.) 

For comparison, we also fit the kerrbb model to simu- 
lated bhspec spectra using the PCA and XRT responses. 
The procedure we follow is identical to that described 
in $3l but we now calculate / co i over the same range 
of a* and / as in Fig. [5l Fig. [6] shows the resulting 
color-temperature corrections for kerrbb fit over a pro- 
portional counter (PCA, solid) and CCD (XRT, dashed) 
bandpass. (We stress again that the different shape of 
kerrbb to that of diskbb means that these are not appli- 
cable to diskbb fits.) Except for the highest luminosities 
of the alpha disk model with a = 0.1, the CCD response 
always yields a higher / co i than the proportional counter. 

A key issue which determines the spectrum is the frac- 
tion of energy dissipated in the region where thermaliza- 
tion will be incomplete i.e. above r e g ~ 3. This crucial 
translucent region of the disk is simply an atmosphere 
(negligible dissipation) for very optically thick disks, as 
the heating occurs at much greater depths. However, 
when the whole disk has r e g < 30 then a substantial 
fraction of the energy can be emitted in the translucent 
region, leading to a large increase in / co i. Fig. [5] shows 
that the alpha disk with a = 0.1 (panel b) is the only 
prescription for which the majority of the flux can be 
produced in the region with r e g < 3 (i.e. logr e ff = 0.5). 
This happens only for luminosities near Eddington and 
is only weakly dependent on spin. A comparison with 
Fig. [5] shows that these models are also the ones with 
highest /coi. The other three stress prescriptions (alpha 
disks with a = 0.01, beta and mean disks with a = 0.1) 
tend to give more modest / co i at these I. This compar- 
ison also identifies r g ~ 30 as the point at which dissi- 
pation in the atmosphere becomes important. At lower 
luminosities (log / < —0.5) T e g is larger and / co i is rather 
insensitive to the stress prescription. 

For these low-to-moderate spins (a* < 0.8) the varia- 
tions of /coi and T e s have a much stronger dependence on 
/ than on a* . These spins are probably the ones most ap- 
propriate for the majority of BHB which have low mass 
companions, as there is not enough mass in the com- 
panion for the accreted material to signifi cantly change 
the a ngular momentum of the black hole (|King fc Kolbl 
1999). Thus they fairly accurately reflect the birth spin 
distribution from massive star collapse, for which the 
best estimates give spin a* < 0.75 — 0.9 ( Gammie ct al. 
2004). However, in the rarer high mass X-ray bi naries 
(which probably also form the ULX population e.g. lKinj 
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Fig. 5. — Effective opticai depth at ffux maximum as a function of spin and luminosity. The contour labels are values of log Teg. The 
panels correspond to alpha disks with a = 0.01 (a) and a = 0.1 (b), mean disks with a = 0.1 (c), and beta disks with a = 0.1 (d). Note 
the change in scale of the contours for panel d. There is less than one order of magnitude change in effective optical depth across the whole 
range, so we decrease the spacing to 0.1 dex as opposed to 1 dex for the other panels. 



2008) there is more mass in the companion star, so ac- 
cretion may change the mass and spin of the black hole 
if the lifetime is long enough. Thus we also investigate 
higher spins, but these are computationally problematic 
for alpha disks with a = 0.1 as they become very effec- 
tively optically thin for a* ~ 1 and I ~ 1. Since this 
stress prescription is ruled out by both spectral and vari- 
ability behavior for the BHB, we take only alpha disks 
with a = 0.01 as a guide to the behavior of the other 
dense disks (mean and beta) as a function of spin for 
0.8 < a* < 0.99. Contour plots for f co \ and r e ff are shown 
in Fig. [7] The contours show a slightly stronger depen- 
dence on a* at these higher values, with r e a decreasing 
and /coi increasing as a* approaches unity. Nonethe- 
less, the overall effect is still relatively small compared 
to the much stronger dependence on I, and the conclu- 
sions drawn from Fig. Q]generalize quite broadly to other 
spins. Thus this predicts that any of the dense disk stress 
prescriptions will give disk spectra at high spin with / cc j 



increasing slightly with I in a way which is not signifi- 
cantly different to that for low spin. 

Thus extreme spins should simply manifest themselves 
as higher temperature disks at a given luminosity, giving 
a clear signature of a maximal Kerr black hole. However, 
none of the high mass X-ray binaries (e.g. Cyg X-l) 
show any signs of this, so we conclude that accretion has 
not yet had time to significantly spin up the black holes 
in these syst ems. The ULXs are also potentially hig h 
spin objects (jEbisawa et all 120031 : iHui fc Krohil 120081 ). 
but these have spectra which are often more complex 
than a simple disk model (see, however, IHui &; Krolikl 
2008), in which case they do not give a straightforward 
diagn ostic of black hole mass and spin (jDone fc Kubotal 
f2006h . 

5. DISCUSSION 
5.1. Stress as a function of radius 
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Fig. 6. — Color correction at flux maximum for kerrbb as a function of spin and luminosity. The solid and dashed contours are 
computed using the RXTE PCA and Swift XRT responses, respectively. See text for more details. The panels correspond to different stress 
prescriptions as described in Fig. [5] 



The observation that stable disk spectra are seen span- 
ning I r~j 0.05 — 0.5 shows that disks are not subject to 
the limit cycles predicted by the radiation pressure insta- 
bility in this range. This rules out an alpha type stress, 
and instead requires that the surface density increases 
(or remains constant in the limit of marginal stability) 
as a function of I at each radius, and this in itself im- 
plies that the accretion disk should remain effectively 
optically thick at high luminosities. Therefore the color- 
temperature corrections will be relatively constant, so 
predicting the approximate L oc T 4 relation for the disk 
dominated spectra. The observation that most binaries 
show an approximate L oc T 4 is strong confirmation of 
this conclusion. 

Sadly, this also rules out the radiation pressure in- 
stability as being the physical mechanism for any ob- 
served behavior for I < 0.5. This is otherwise an attrac- 
tive possibility for the origin of the very high state, the 
other spectral state seen in binaries at high mass accre- 
tion rates. Although this does not show the predicted 



limit cycles, there is the possibility that the non-linear 
outcome of the instability could in stead lead to strong 
Comptonization of the disk (e.g . [Rozarisk a fe Czernvl 
l2000HKubota fc Makishimall2004f) . Nonetheless, the ar- 
guments above preclude an alpha type stress, at least 
below I ~ 0.5, yet the very high state is indeed seen in 
the s ame range in luminosity as the disk dominated state 

(e.g. lRemillard fe McClinto^l200^ 1. 

At higher luminosities (up to I ~ 3. lDone et al.l f2004). 
there is a limit cycle observed in GRS1915+104 which 
could be the remnant of this radiation pressure instabil- 
ity from a mar ginally stable stress such as predicted with 
the mean disk (jHonma et"aLlll991b iMerloni fc Navakshinl 
l2006f) . However, this limit cycle behavior could alterna- 
tively be from other instabilities to do with the super 
Eddington flows. This is an important distinction, as 
the classic radiation pressure instability is not scale in- 
variant with mass, but is triggered at / oc to -1 / 8 . Thus 
even if the disk is stable in stellar mass black holes up 
to I ~ 0.5, it can be unstable at I > 0.09 for a 10 7 M Q 




Fig. 7. — As in Fig. \E\ but extended to high spin for an alpha disk with a = 0.01. The left panel gives effective optical depth at 
flux maximum, while the right panel gives the values of f co \ for kerrbb computed using the RXTE PCA (solid) and Swift XRT (dashed) 
responses, respectively. 



AGN. In this case the instability could play a role in 
producing the puzzling, 'soft excess' seen in high mass 
acc retion rate AGN, predominantly Narrow Line Seyfert 
l's (jRozariska fc Czerny||20Q0t iFabian et al.ll2002ft . How- 
ever, there are multiple similarities between the very high 
state and Narrow Line Seyfert 1 spectra to make it more 
likely that there is a similar explanation for both types 
of object. 

Similarly, those properties of super Eddington flows 
which likewise depend on the alpha stress are also prob- 
ably not a realistic description of these disks. Radiation 
trapping is a generic feature of any of the stress pre- 
scriptions fe.g. lHonma et alii 199 If ), resulting in optically 
thick advection of en ergy especially in the disk midplanc 
(|Watarai et al.lr2001h . However, the denser disks which 
result from the alternative stress prescriptions are less 
likely to become effectively optically thin to the escaping 
radiation, so are unlikely to show the very large color- 
temperature correc tion which can arise from overheatin g 
of alpha disks fe.g. lBeloborodovl ll998; K awaguchll2003f ) . 
We caution that fitting such models to high mass ac- 
cretion rate spectra from ULX and NLSl's may not be 
appropriate. 

5.2. Caveats 

So what then are the caveats to using disk spectra as 
an estimator for black hole spin? The first is that the 
spectra must be disk dominated. The disk radius can be 
under - or overestimated when a substantial fraction of 
the dissip ation goes instead into a corona, either in the 
very high (jKubota fe Dondl2004t iDone fe Kubotall2006f) 
or low/hard state (GDP08). This may be the origin of 
the discrepancy in spin determina tion by disk spectral fit- 
ting in GRS 1915+104 (compare iMiddleton et all 120061 : 
iMcClintock eTaII[2006l) . The spectra of GRS 1915+104 
may well be more complex than the disk dominated 
spectra seen from BHB due to its higher luminosity 
(|Done et al.ll2004D . perhaps powering strong winds (see 
below). It seems premature to apply the uncertainties 



from this one pathological object to other sub-Eddington 
BHB fe.g. lRevnolds fe Fabianl[200cl . 

Secondly, the value of the color-temperature correction 
is robust for all of the dense disk stress prescriptions, but 
varies with bandpass and also varies with the detailed 
disk model used to fit the spectra. We caution that dif- 
ferences in spectral shape between the models means that 
the same color-temperature correction factors cannot be 
simply applied to diskbb and kerrbb fits. Ironically, for 
both the generic proportional counter and CCD band- 
passes, the diskbb fits give a simpler (approximately 
constant color temperature) representation of the data. 

Thirdly, the structure of the disk as I — > 1 is not well 
modeled by this code. We neglect optically thick ad- 
vection, and changes in the scale height of the disk with 
radius which may lead to self shie lding for highly inclined 
objects (e.g. IWatarai et aLl l2005), but these are unlikely 
to have much effect for I < 0.5, where the majority of 
BHB data are taken. 

However, disks at high luminosities are likely to also 
power substantial winds. There are clear signs of an 
equatorial disk wind from observations of highly ionized, 
blucshiftcd, Fe Ka absorption lines. These are ubiqui- 
tously seen in CCD data from high/soft state BHB (i.e. 
I > 0.1) with inclination angles > 60° (e.g. DGK07). 
This material is probably a thermally driven disk wind, 
where irradiation of the outer disk layers produces a 
Compton heated skin whose high temperature allows gas 
particles to have enou g h energy to escape at large radii 
(|Begelman et al.l Il983t IWoods et al.l 119961) , though this 
may be (substant ially?) enhanced by magnetic fields 
(|Miller et a l. 2006) . Such winds almost certainly become 
more powerful with I, especially as / — > 1 and radiation 
drivin g enhances them still further (|Proga fc Kallmanl 
2002). Counterintuitively, these stronger winds may be- 
come less observable as the stronger irradiation means 
that they become more ioniz ed, losing even the He- and 
H-like Fe (|Woods et al.lll996h . Thus there may be winds 
which have substantial column, but are unobservable as 
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they are completely ionized. Electron scattering has the 
effect of suppressing the observed luminosity by a fac- 
tor e _T for lines of sight intersecting the wind, while the 
scattered flux (1 — e~ r ) can enhance the luminosity along 
the disk axis by a factor 1 + (1 — e~ T )Q,/Aii where f2 is 
the solid angle subtended by the wind. This may be the 
origin of the occasionally observed 'bends' in the L — T 
diagram, where the data deviate away from the generally 
observed slow increase in fco\ as a function of luminosity 
e.g. in GRO J1655-40 there is a much stronger increase in 
/col with I than predicted by the dense disk models (e.g. 
Fig lb), while in LMC X-3 there is a slight decrease (see 
e.g. lDavis et~alll2006| ). While the inclination of LMC X-3 
is not well known, GRO J 1655-40 is plainly at very high 
inclination, and ha s clear wind features detected in m ul- 
tiple observations (jMiller et al.ll2006fc ISala , et alJl2007f ). 

Thus as well as avoiding spectra with strong tails, and 
super Eddington sources, we also caution against any 
serious attempt to derive spin from 'bent' L—T diagrams, 
where winds may be important. 

Our results also hinge on the fact that the dissipation 
above the effective photosphere is small when the surface 
density is sufficiently large, so our assumptions about 
the dissipation profile wit h height may be an important 
cavea t. Previous work (|Davis et al.l [2005; Blac s etaLl 
120061) has demonstrated that dissipation profiles moti- 
yated by rad i ative, stratified sh earing box simulations 
(|Turnerll200l iHirose et al"1l2006l ) do not put enough dis- 
sipation near the surface to significantly alter the spec- 
trum, although a more exhaustive range of parameter 
space needs to be explored. However, the alternative 
stress prescriptions discussed here are partially moti- 
vated by theoretical arguments that limit the magnetic 
stre sses due to transport of buoyant ma gnetic field (see 
e.g. lStella h Rosner|[l98l lMerlonill2003D . If transport of 
magnetic field is significantly larger than currently seen 
in simulations, dissipation would be concentrat ed much 
neare r the surface than our models assume (|Merlonil 
2003). This would give significant dissipation above the 
photosphere, and hence lead to a marked increase in 
color-temperature correction compared to those calcu- 
lated here for the dense disks. As long as this is not a 
strong function of luminosity then it could still fit the 
data, but would lead to the spin being overestimated. 
Since the spins are generally low-to-moderate, there is 
not much scope for the real stresses to produce a signifi- 
cant fraction of the dissipation within the photosphere. 

Similarly, the data also preclude a large amount of dis- 
sipation at the last stable orbit. If this were to thermalize 
it would lead to higher temperature emission, so current 
models would overestimate the spin, yet these already 
favor low spins. Alternatively, if the dissipation did not 
thermalize, it would instead produce a substantial non- 
thermal component in the spectrum, yet the spectra are 
dominated by the disk emission in the high/soft state. 
Thus there is only limited scope for stresses at the last 
stable orbit in disk dominated spectra, in contrast to the 
results from non-radiati ve, relativistic magnet ohydrody- 
namic simulations (e.g. iBeckwith et al.ll2008l ). It may 
be that the character of the flow qualitatively changes 
when radiative cooling is included in the simulations, as 
the small scale height fields sustained by a thin (radia- 
tively cooled) disk may lead to much less stress at the 
last stable orbit than the large scale fields which can be 



generated in the large scale height, non-radiat ive flows 
(jAfshordi fc Paczvhskill2003l : iShafee et al.ll2008l K 

This potential difference in dissipation between large 
scale height flows and thin disks may actually be ob- 
served in the behavior of the transition between the 
low/hard state, where the spectrum is produced by 
Compton scattering in a hot, optically flow, and the disk 
dominated state where the emission is quasi-thermal. 
The low/hard state is probably produced by a large scale 
height flow, which has an efficiency which is at least a 
factor 3 lower than that of a (stress-free inner bound- 
ary condition) thin disk down to the same radius (e.g. 
iNarayan &: Yil[T995h . This predicts that the luminosity 
at the transition should jump by at least a factor 3, yet 
the data show a smooth transition of less than a factor 
~ 2 (e.g. th e compilation black hole binaries, especially 
Cyg X-l, in iDone h Gierlihskil [200l . This can be ex- 
plained if the large scale height, less efficient flow in the 
low/hard state has continuous dissipation across the last 
stable orbit due to large scale height magnetic fields. In 
effect this allows the flow to extend down to smaller radii, 
so it taps a larger fraction of the gravitational potential 
energy, hence somewhat compensating for its lower radia- 
tive efficiency compared to that of the thin disk (Done 
et al. 2008, in preparation). 

6. CONCLUSIONS 

The detailed nature of angular momentum transport 
remains a significant uncertainty in our understanding 
of accretion flows and their emission. It is not yet en- 
tirely certain that the MRI, which is the best candi- 
date for the source of the turbulent stress, truly does 
saturate at a level high enough to satisfy observational 
const r aints (iPessah et al.l 120071 : iFromang &: Papaloizoul 
120071: iKing et alj I2007D . Even if this question is re- 
solved in favor of the MRI, a detailed understanding 
of the stress and resulting disk structure will likely re- 
quire fully-relativistic, global simulations with realistic 
thermodynamics. Since such simulations require signifi- 
cant advances in computing power, we are forced rely on 
much simpler parameterizations of the stress until these 
become available. Even though these prescriptions are 
ad hoc, they may still capture important aspects of the 
physics in real astrophysical flows and provide useful con- 
straints by direct comparison with observations. 

With this motivation, we have computed the accretion 
disk spectra predicted by several different stress prescrip- 
tions using the most complete spectral code currently 
available. This includes both non-LTE ion populations, 
with radiative transfer in the disk and full general rel- 
ativistic ray tracing to propagate this flux to the ob- 
server. These models give results which are very close 
to the observed behavior of the spectra of black hole bi- 
naries in the thermal dominant (high/soft) state from 
L/L,Edd ~ 0.06 — > 0.6. The majority of this data comes 
from proportional counters, and generally shows that the 
temperature changes slightly more rapidly with luminos- 
ity than expected for a disk of constant inner radius, and 
constant color-temperature correction (e.g. the compila- 
tion of DG04 and DGK07). This is as predicted by all 
the high surface density disk models presented here, and 
should hold generally for all stress prescriptions where 
less than 10 per cent of the energy is dissipated above 
the effective photosphere. The surface layers then sim- 
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ply act like a passive atmosphere, with properties set by 
the effective temperature. The increase in temperature 
with luminosity gives rise to a small increase in color- 
temperature correction and the resulting spectra are re- 
markably similar irrespective of the detailed form of the 
stress. 

Such model fitting can also place interesting con- 
straints on the angular momentum transport. The clas- 
sic alpha disk with a = 0.1 becomes effectively optically 
thin at the highest luminosities. The resulting increase 
in color temperature with luminosity is much more rapid 
than observed from existing data, ruling out such stress 
prescriptions. The standard alpha disk is also unstable to 
the thermal-viscous radiation pressure instability which 
predicts limit cycle behavior which is not observed, again 
showing that the disk is denser than such models predict. 
However, even these models give the same spectra at low 
luminosities (L < O.lZ/Edd), where the surface density is 
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sufficiently large for the disk to remain very effectively 
optically thick. Therefore, where we have a clear view of 
the disk, unaffected by a moderately optically thick disk 
wind, the disk spectra should provide a relatively robust 
estimator of the disk inner radius, and plausibly, the spin 
of the black hole. 
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